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ABSTRACT 


A possible alternative to the ordinary gas-like computation 
for nuclear matter is investigated under the assumption that the 
nucleons are arranged in a lattice. BCC, FCC and HCP structures 
are investigated. Only HCP shows a minimum in the energy vs. 
density curve with a modest binding energy of - 1. 5 MeV. The very 
low density limit is investigated and sensible results are. obtained 
only if the tensor force decreases with the density. A study of the 
elastic properties indicates that the previous structures are 
mechanically unstable against shearing stresses. 



I. INTRODUCTION 


It can hardly be denied that the problem of nuclear matter 
represents one of the most difficult and at the same time the most 
challenging problem in nuclear theory. The principal reason for 
rather an unsatisfactory solution of this central problem in nuclear 
physics is that, one's attempts to reproduce the value of the volume 
energy as close to 17 MeV as possible, have been frustrated largely 
due to successive improvement of the theory with the tendency now- 
adays to stabilize around 11. 1 MeV. This value was computed by 
Siemens^ using the most sophisticated many- body theory now 
available, namely, the BBP approach. A long way (about 6 MeV) 

remains to be accounted for, however. In a recent thorough review 

(3) 

of the subject, H. Bethe has indicated how contributions of a 
different nature (none of them being absolutely watertight) can be 
added to account for just about the missing 5 “ 6 MeV. Apart from 
the uneasiness felt by the process of adding to a good solid number like 
XI. l MeV, a less firmly established contribution, the crucial question 
remains as to whether the theory is actually reliable enough to be 
considered the appropriate basis for understanding the shell model. 
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The 5—6 missing MeV are too big a chunk to be taken lightly. On 
the other hand, one could assume the following dichotomic posture: 
(a) the sought for volume energy is not actually 17 MeV but much 
less, say about 12 MeV. This possibility was proposed some years 
ago by one of the authors' by showing that summing the energy of 
a finite nucleus, 

E <«|T + |v|a> = iz e a + l £ <a l T l a> 

a a a 

upon using a reasonable set of single particle energy levels, an 

analytic formula was obtained which could in turn be rewritten to 

exhibit the dependence of the mass formula on A and Z. All the 

well known terms, namely, volume, surface, asymmetry and 

Coulomb were obtained. A new term was, however, found to 

appear with the same sign as the volume energy, whose dependence 

-2/3 

turned out to be A’ . The term was shown to contribute 3^4 
MeV for medium heavy nuclei. The various coefficients were given 
analytically as functions of parameters usually encountered in the 
shell model or optical model computations. The overall signifi- 
cance of the computation does not rely on its numerical aspect, 
but rather on showing that there could still be room for discussion 
about the, proper limit of E/A where A -* » . (b) A second line of 



5 


( 5 ) 

thought was first introduced by Overhauser. He pointed out that 
a common feature of any many- body treatment is the yet unproven 
hypothesis that the ground state of a Fermi gas is the familiar 
sphere of occupied momentum states, instead of, for instance, 
spatially periodic Hartree-Fock solutions, It therefore remains to 
be seen if any other type of a single particle wave function could in 
principle produce a lower energy. Plane wave functions bring out 
the flavor of a gas, whereas, the abstract concept of nuclear matter 
could also be thought of as a solid or a liquid. This paper reports 
the results of a full many- body computation intended to indicate 
that the nuclear matter is indeed not a solid or a quantum crystal. 

If it is possible to make a definitive statement about nuclear matter 
not being in a crystalline state, any extrapolation to liquid structure 
is of course fraught with many more uncertainties. Liquids are 
in general, difficult to distinguish from solids in many respects. 
Liquid Ar, for instance, clearly shows a crystalling structure 
in its first shell with about 8 neighbors instead of the usual 12 in 
the close packed structure, when viewed with neutron scattering 
experiments. Other liquids still show solid- like structure with an 
increased number of nearest neighbors. The whole concept of a 
liquid is not clearly under stood, as indicated by the partial success 
achieved by distribution function theories and by Lennard- Jones- 
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Devonshire theory that reproduce reasonably well the data, even 
though their starting point is quite different. The only clear-cut 
definition upon which nobody disagrees is that liquids do not 
possess any ability to resist shear stresses whereas solids do. 

This mechanical definition is of fundamental importance in 
the following sense. Suppose one computes the ground state energy 
of a many- body system of nucleons by using, instead of plane waves, 
wave functions such as Bloch wave functions or Wannier functions 
and finds that the energy at the nuclear density turns out to be less 
than the one yielded by the corresponding gas-like computation. 
From that result it would be quite premature to conclude that the 
solid structure necessarily represents the actual state of affairs. 

The mechanical properties of the solid have to be tested. The 
lower energy state can turn out to be totally unstable against 

f 

shearing stresses, i. e. , one has to make sure that the elastic 

( 8 ) 

constants of the medium do not go soft. M. Born has given a 
series of criteria to decide as to whether the ordered structure 
melts, evaporates or goes into an intermediate state similar to 
plastic or "gel. 11 

We have attempted three different ordered structures: BCC 

(body- centered cubic), FCC {face- centered cubic) and HCP 
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(hexagonal close packing). In the three cases, the system was 

shown to be mechanically unstable against shearing stresses, or 

put in a different way, the significant elastic constant went soft 

12 14 

for the density integral between 10 gr/cc ^ p ^4x10 gr/cc. 

One can therefore conclude quite safely that the idealized N = Z 

14 3 

infinite system at p =s 3. 10 gm/cm is surely not in any of the 

most commonly encountered solid state structure. On energetic 

grounds it turns out that the binding energy per particle E/N has 

14 

a minimum of - 1. 5 MeV at 2. 842 10 gr/ cc for the HCP structure 
whereas for BCC and FCC the problem of the existence of a 
minimum is much harder to discern. 


II. THE MANY- BODY TREATMENT 

The many- body treatment used is an extension of the T-matrix 

3 

approach to quantum crystals (solid He ) recently discussed by 
Brandow. ^ An excellent review can be found in a paper by Guyer^^ 
The Slater determinant for a system of N particles is built up by 
single particle wave functions of the form 

3/2 -3/4 -(a 2 /2)|r-R| Z 

0{r) = a tt e /it ( 1) 

where R is the coordinate of the lattice site around which the 


particle is supposed to perform an oscillatory motion under the 
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combined influence of the remaining (N- 1) particles. Equation (1) 
is the eigenfunction of an harmonic oscialltor potential U(r), 
centered around the lattice 3ite R, i, e. , 

U(r) = — muu 2 (r-R) 2 . (2' 

The frequency <d (which cannot exceed the Debye frequency oj^ at 

zero temperature) enters in the wave function (1) through the para- 
- 1 9 

meter a = (irwM) s that represents the spread of 0 around the lattice 
site. At the present stage of the computation such a parameter is 
as yet unknown and plays a fundamental role in the theory. Some 
remarks about the explicit form (1) are quite in order. One is 
accustomed to think that a wave function of a particle in a crystal 
should possess the periodicity of such a crystal or equivalently be 
a Bloch function. To begin with, the HF equations are highly non- 
linear and if one starts with a 0 of that form one will not necessarily 
end up with such a form after the desired convergence has been 
achieved. In the second place, any linear combination of Bloch 
functions is still a solution; such a linear combination can turn out 
to possess quite a high degree of localization (such as Wannier func- 
tions) of the form given in eq. (1). Evidently one of the prices paid 
by not actually performing such a linear combination of Bloch functions 
is represented by the fact that the parameter u) is unknown. The 
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troubles encountered in not making such a consideration and bodily 

considering 0 as a trial function with a as a variational parameter 

3 

are too well known in solid He not to be repeated here. Nosanow 

( 12 ) 

and co-worker have determined a from HF. We have followed 

the same spirit in our calculation, i. e. , we take the single particle 
potential to be given by 


U(r x ) = ZJtfj < r 2 > V (r r r 2> **j (r 2' d r 2 


where the two- body nucleon- nucleon potential is taken to be Reid's 
phenomenological soft-core potential. The sum over j indicates 
the sum over the successive neighbor distances (shells) from a 
given particle located at R-^ in an arbitrary system of coordinates. 
The values of ijo (or a) obtained by solving eq. (3) are then used to 
start the full HF equations which finally give the final form of 0. 
The angular momentum decomposition of Vtr^-r^) complicates life 
quite a bit and it was found that. the following decomposition was 
very useful 


imR 
e 

R 


Jfc £ < 2 ^ 1 > J t + l/2 (mp)H ^ + l/2 (mr)P i <COS CP) 


p < r 
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Reid's potential is in its most important part a superposition of 
Yukawa's. The structural nature of the specific crystal BCC, FCC 
or HCP enters when the summation over j has to be performed. 
Depending upon the density, the sum was performed over as many 
neighbors as necessary to ensure a satisfactory convergence. At 
nuclear density the number is about 5—6, whereas at much higher 
densities, say p ^ 10 gr/cc, one has to sum over at least 36 
neighbors. The values of a , the spread of the single particle 
wave function vs. p (density) are given in Tables 1, 3, 5 together 
with the first neighbor distance. Each particle spreads 
about 30% of the distance of the lattice site, an excursion quite 
acceptable even on the basis of an empirical rule such as the one 
often encountered: The crystal melts whenever the oscillations 

reaches a sizable fraction of the lattice constant (Lindeman rule). 
Even though this criterion can be easily seen not to have general 
validity (to have a semi- quantitive basis, zero point oscilllations 
are very large and yet do not induce any melting)our spread is well 
within the safety margin for each nucleon to perform comfortable 
excursion from its lattice site. We should perhaps remark that 
Lindeman’ s rule can be derived in many alternative ways. 
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A major complication arises in the nuclear case because of the 
possible spin configurations: given a structure, say BCC, one could 
in principle arrange the protons and neutrons in pr actically an 
infinite number of ways depending upon the spin configuration. In 
a BCC structure one has two interpenetrating cubes amounting to 
16 site s t 8 for neutrons and 8 for protons; the configuration is still 
infinitely degenerate with respect to spin arrangements. One could 
start with any site of a cube with, say, a neutron with its spinip; 

in the next site on the same cube one can put Nt , Nl, Pi or PI. In 

the next point one again has the same degeneracy and so on. One 
can come up with the most asymmetrical configuration with three 
consecutive particles with spin t, then two with spin 1 and thus 
randomly distributed as far as the unit cell is concerned. We have 
tried many of them with an arbitrary degree of randomness just to 
insure that they are not energetically convenient. It was then estab- 
lished that this was indeed the case and that: for a given configuration 
the minimum was achieved whenever the spins (of neighboring parti- 
cles on the same cube) were symmetrically arranged, half with spin 

ou- 

T* 

up and half with spin down. 


Natura non facit saltus (Lucretius, de Rerum Natura). 
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An analagous situation was found to hold for the other two 
configurations studied, namely FCC and HCP. In the usual nuclear 
matter computation one sums over all the possible spin orientations, 
an operation that is not allowed in a discretized structure such as 

an ordered arrangement. This will explain why the energy per 
particle formula looks different from the one that is usually 
displayed in Brueckner formulation. 

The many- body formulation used here is the same that has 

3 

been quite successfully applied to solid He , or for that matter, to any 
quantum crystal. The main difference will be in the angular momen- 
tum decomposition of the two- body wave' function, say tjf. In the 
3 

solid He case, the formulation includes only 1=0 waves since the 

Leonard- Jones potential is assumed to be a 1=0 state. This is not 

( 12 ) 

at all obvious and the authors have recently shown that a much 

better agreement with the experimental data can be achieved if one 
3 

allows the He atoms to interact also in 1=1, 2 states. In principle, 
it is difficult to fully understand a wave function decomposed into 
angular momentum waves upon considering that the Y m are related 
to a sphere, whereas the particles are located spatially with a cubic 
symmetry for BCC and FCC, and hexagonal symmetry for HCP. On 
the other hand the phenomenological nucleon- nucleon potential now 
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2S+ 1 

in extensive use is given explicitly for each L wave and one is 

J 

therefore forced to decompose the wave function into Y 1 * 1 . About 

'L 

this point one could sophisticate a little bit more by considering how 
the potential itself is deduced from the S-matrix. Reid potential is 
not strictly under this jurisdiction since it is supposed to have been 
deduced upon fitting scattering data and deuteron properties. Even then 
the general basic- structure is written down on the pionic theory. When 


one deduces V(x), one essentially Fourier- transforms the second 
order of the S-matrix which in turn makes use of the Dirac free 
spinors. Evidently this is not the case in a crystal structure where 
the nucleons would have a relativistic harmonic- oscillator wave 
function (which do not exist) as a spinor. The use of a free two- 
nucleon potential could in principle be questionable but this is the 
best one can do at the present time. 


The actual computation of E/N has been reviewed elsewhere 
and we will therefore only quote the results. Since the original 
Slater determinant is made up of displaced harmonic oscillator 
wave functions, the kinetic energy will evidently be 


E 


3 ft' 


K. E. 


m 


2 

a N 


( 4 ) 
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The potential energy has a familiar form too, namely, 


3 j 3 

d x. 


E 


P, E. 


^ J$(ij) V (ij) Y (ij) d Xj 

3 3 

i<j (ij) Y (ij) d Xj d x. 


(5) 


where $ is the free two- body wave function, i. e. , $(r^r^) - <b{r , 

whereas Y(r r ) is the correlated or perturbed wave function which 
1 2 


we have to study. Transforming to relative and center of mass coordi- 
nates, expression (5) becomes 


J F. E. 


= _L nS f»(r)V(r) T ^ L dl r = 1 nE 


(r) f (r) d r 


( 6 ) 


where N is the number of particles at distance A from the one 
A 

chosen as the origin. In any non-T-matrix approach, like in 

the Jastrow variational method, one does not have at this point an 

equation for f of the same form as the one we are about to present. 

The form of E / N is also different since a logarithmic term is added 

E E. 

to the potential. Among the various degrees of sophistication that 

one can use to write down the equation for ijf , we shall make use of 

(13) 

the one first employed by Guyer and Zane, 1 . e. , 
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T 1 + T 2 + U(l) + U(2) + V (12) 


<jr ( i2) = ef{i2) 


(7) 


which is physically easy to understand. T(l) and T{2) are, 
respectively, the kinetic energies; U(l) and U(2), the harmonic 
oscillator potentials of particles 1 and 2; and V(12), the two- body 
interaction potential. Equation (7^ is the equation of motion of two 
particles, each of them moving in an harmonic potential centered 
around two different points, and, in addition, interacting through 
a two- body potential. When the two harmonic potentials are 
centered at the same point (in par ticular, R 1 = R 2 “ tlie P roblem 
simplifies greatly and it was studied in connection with shell model 
computations by Nigam some years ago. Equation (7) is the 
most difficult equation to solve. Once is known, the energies 
can easily be calculated. If one uses for U(l) and U(2) the 
form (2), equation (7) becomes 


( 8 ) 


{. AL 7 2 . iL v 2 + i u) 2 / m r^ + m r 2 \ + \ uu-V - uu> 2 ?A + V(r){f = M 
2M R 2p r 2 \ 1 1 2 2/2 l 


where jjl = 


m i m 2 

m l+m 2 


, M = m 1 + m^ , r " r 1 ’ r 2 * R = 


m i r l + r 2 m 2 

M . 


At this point it is usually assumed that Y = 0(H0(2)g(r): if, 

on one side, it is quite a pictorial way of stating that Y differs from 


$ by a g(r) (correlation function) , on the other side, it has quite 

a disastrous consequence if one filters $(1)0(2) to the left side of 

eq. (8) (as is usually done), leaving an equation only for g(r). 

Equation (8) as such, contains, after introducing Y* = Y/r, only 

second derivatives, whereas the process of filtering 0{1)0(2) to 

the left, introduces a first derivative which cannot be eliminated: 

from the numerical point of view this is highly undesirable, as 

we have found by personal experience. As a matter of fact, the 

only computation that has made use of the equation for g(rl instead 
( 15) 

of for Y itself, shows that the boundary conditions were more 

difficult to meet. If, moreover, the potential V(r) has an 
angular dependence, then the transportation of 0(1) 0(2) to the 
left is totally illegal, since the remaining g(r) would still have to 
be decomposed into partial waves, thus, 


g(r) = 2 Y ^(n) 


( 9 ) 


without the possibility of giving any physical significance to t 
whatsoever, since part of the physical angular momentum (most of it!) 
contained in 0(1) 0(2) has been transported to the left. To be more 
specific, a matrix element of the form 
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with t and £ defined by eq. (9) , would be completely meaningless. 
Since in the nuclear case, V(r) is highly angular- momentum dependent, 

the only choice is to write the full angular momentum decomposition 
of' t as 


CO I , 

^ = ,?o ^ (U+ V (r ' ^ 


( 10 ) 

in obvious notation. The second summation over is introduced 
to take into account the presence of tensor forces. Equation (8) 
contains a cos 0 term which (much like the Stark effect* couples l 
with l ± 1. It does not couple different spins and therefore it is 
to be expected that S^, will be correlated by such a term, 

and analogously the triplet states 3 S 3 p 

1 1 O’ 1, 2’ 

3 3 ~ 3 

^ 1 ’ ^2, de linked together through the solid state term. 

After substituting eq. (10) into eq.-(8) and performing a long angular 


momentum algebra, the results are (x = r/r , d = A/r 


O’ 


a = ar Q ): 


SINGLET 


\ + (E - U 0 )h Q -F a(x,d)h 1 


+ (E - u 1 ^ h 1 + 3’ a ( x * d ) [h Q - Zh 2 ] 


h" + (E- U 2 )h z - — a(x, d)h 


= 0 , 


( 11 ) 
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where the following notation has been used 


00 , 

r ho f h o 


, 10 u 

hi =h i 


20 , 
r + 22 = h 2 


a(x, d) = — a^xd 


TRIPLET 


__ 2^ r . . 

h" + (E-U,)h - 78 ^ h, + ± a(x, d) h,,+ 3h„+5h, 

5 oi T . 2 3 9 456 

ft 


~ ~ _ 2ur 

h 3 + < E - U 3 )h 3 ' ^ V T — h 3 ' 9 a(X ’ d> h 4"i h 5 + 2 h 6 

ft 


h"+ (E-U 4 )h 4 + 


h^+ (E-U 5 )h 5 + 


h 6 + < E - U 6 )h 6 + 


0 +I 


+ ^ a(x, d) (h 3 + 72 h ? j- 72(h 3 + 72 =0 

~7 a(x t d) |h + hj- ( h + -—;h \ - 
1\ 5a/ 2 / 5a/2 \ 5a/2 7 


9 2 

h 7 + < E - U 7 )h 7 - 7/8 V T-T h 7 ' 9 a(Xid) 2h 4 + K + iK 

* L J 


h"+ (E-U 7 )h ? - 78 V T — 9 


2 Jf0 ■ f^a(x,d) h 4 -|h 5+ ih 6 


0 - — a(x,d) 3h c + h, 

10 5 6 


( 12 ) 
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with the following notation 



11 

r ^oo 

= h 3 

3 <v 

21 I. 

r hi = h 6 


3 Sr 

11 

r *02 

= h 3 

3 °i : 

11 u 

r\!/ = h_ 


3 *v 

01 

r *ll 

= h 4 

% 

ll r 
20 = h 7 


3p r 

11 

r hi 

h 5 


21 . 
^22= h 8 

(13) 



As one can see, the equations are intrinsically coupled by 
the solid state term a(x,d'. Only two normalization conditions are 
therefore necessary: one for the singlets and one for the triplets. 
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i 

The solution of the 3+8 equations has presented us with major 
numerical difficulties; no simple theorem of the Sturm- Li ouvi lie 
form exist for coupled equations and the search for the ground state 
energy was indeed a painful and delicate process. 

The major concern was to make sure that the h^/s ^ad no 
nodes or else that the energy eigenvalues being computed were actually 
the lowest. In problems of stellar evolution and in studies of turbu-. 

fence, equations of this nature frequently arise, but with the major 
difference that seldom any eigenvalue problem has to be handled; 
the energy E is just a prescribed parameter. 

-1 

Given a certain A (first neighbor distance) and a (spread of 

t 

the wave function) we found that around the nuclear density 


t 

p = nm, 


4 2 fZ # 

3 3 ’ 3 * r FCC “ r 0 

r FCC r BCC r HCF 


the 3+8 solutions were behaving very well; by that we mean that 
the maximum of the h' s was at the proper neighboring distance, A . 
The shape of h^_ is essentially a gaussian as seen in the figures 
below; this is to be expected since Y differs from $ by a correlation 
function that doesn't usually change the shape. The location of the 
h f s is very important since the energy of the system is proportional 
to hV; classically hV is just V computed at exactly the lattice 
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site since the wave functions are all delta functions, i. e. , when 
A -> 00 > <t> -» ^ {fj-Rj). Quantum mechanically, besides the location 
of the peak determined by the 3+8 equation, each h has a certain 
spread, a . For the same r , if A increases the h’s have a 
tendency to lag behind. 

We have tried three possible configurations: BCC, FCC and 
HCP. The geometrical disposition of neutrons and protons as well 
as the energy per particle are given below for each case. 


Energy Expressions 


As mentioned earlier, the energy expressions depend on the 
spin arrangements. We will therefore have to distinguish among 
several cases. Using Eq. (6) and the definition of \|r as given by 
(10), the various possibilities are 


! 

i 


i 
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« N +# N+ (OR pt pt) 

{ a(x)j l{ h 6 V 6* h 5 V 5 } dx 

e l = 

' Ja(x)j 1 [h 6 + h 5 ] dx 

Nt/ N + (OR pt P+) 

0> 

' : r - J a <*>£ j 0 V O h O + 5j 2 V 2 h zl' k f a(x »l! 2h 6 V 6 + h 4 V 4 ] dx 

h = A ~ +o~a>— 

* . • J a ( x )[j 0 h 0 +5j 2 h 2 ]dx | afxjj^^+h^dx 
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As explained before, only two normalization constants are required. 
They are simply determined by requiring that each energy denominator 
should be normalized to unity in conformity with the man- body theory 
requirement that 

_[$. fd 3 r = 1 . 

The explicit forms of the constants A, B, C, D, E, F are finally obtained 
and since they do not present any interest per se they have not been 
written down here. 


The remaining quantities are defined as 


a(x) = xe 


-a x /4 ; 


1 2 

y = — a xd 


j 0~y Sinhy ; coshy - -y sinhy ; j 2 = -~coshy - (“ + s i 

V V ' 


sinhy 


The origin of the Bessel functions (of imaginary argument) can be 
easily understood. When, in Eq. (5), we write the product of 
0(1 )0(2) for ${ij) the result contains a term of the form 


exp 


( 4 ?- &) 


which, when expanded in plane wave in a form analogous to Eq. (10), 
gives 


X 


M S e (a 2 /2)?.A = £ i X V4n ( 2Hl> (?• a) 

\ f J 


24 


if only 4 0 one easily recovers the expression for the energy 

( 16 ) . 3 

used by Massey and Woo in their study of solid He . 

If we sum up all the energy expressions, i. e. , we sum 
over the spin orientation, we recover the familiar expression of 
Brueckner theory. 


III. RESULTS 
B. C.C. Structure 

The geometrical disposition of the nucleons is shown in 
Figure 1. Out of a plethora of possibilities, the one given 

below was found to minimize the energy. In Table 1 we present 
the results: the first column gives the density; the second 

the appropriate lattice parameter; the third the location of the 
first neighbor; the fourth the spread of the wave-functions; then 
the kinetic, potential and total energy per particle follow in order. 



TABLE 1 


p xio 

r BCC 

A 

-1 

a 

KE 

PE 

E/N 

(g/cc) 

<f) 

(f) 

<£) 

(MeV) 

(MeV) 

(MeV) 

1. 045 

3. 175 

2. 75 

1. 04 

28. 63 

-28. 90 

-0. 2744 

2. 842 

2. 274 

1. 97 

0. 73 

57. 61 

-54. 95 

+ 2. 667 

3. 96 

2. 037 

1.76 

0. 65 

74. 18 

-66. 54 

4 7. 642 


14 3 

At nuclear density p = 2. 842 x 10 gr/cm , the system 
is unbound with E/N = + 2. 667 MeV. To appreciate the nature 
of the wave- functions we. list below in Table 2 the distance 
in fermis where they peak at the various densities and the 
first neighbor distance around which they should peak. 




O' 


TABLE 2 

LOCATION OF THE PEAK OF THE WAVE FUNCTIONS (in fermis) 


p X 10 14 
igf cc) 

r _ _ 
BCC 

it) 

ls o 

V 

1 


A 

(£) 

3 s 

S 1 

3~ 3 

S P 

1 0 

3 3 

P 1 L 

3 3~ 3 

1 Dl 2 

1.045 

3. 175 

2. 28 

2. 79 

2. 81 

2. 74 

1. 651 

1.64 



3. 96 

2. 037 

1.63 

1. 87 

1. 95 

1. 76 

1. 38 

1.14 1.79 

1.71 1.73 

1.4 1.42 1.72 
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preceding page blank not filmed 


Here, too, at the nuclear density p - 2. 84 x 10 gr/cm , the 

system is unbound with an energy of + . 5936 MeV, We have also 

added the values of E/N for much lower p where the peak h* s occur 

at an nonphysical distance. It is seen that at r = 10f, or p = 

FCG 

12 3 

6.69 x 10 gr/cm , the potential energy is -8,01 which is clearly 
absurd,. The reason can be easily understood by looking at Table 4, 
where we present the location of the peaks of the various h r s 
(compared to where it should) for different densities. The situa- 
tion becomes disastrous at lower densities and this explains the 
unrealistic negative energy. 

In Figures 3—18 we display the h' s at r p££ ~ 2. 566 for 
two distances A = 1.8 If (first shell), and A = 9. 07f (24th shell), 

and also exhibit the corresponding potentials in the combination 
hV. It is seen that the location of the peak is excellent for 
A = 1.8 If, whereas at A =9.07, the h’ s lag behind. This is 
of no particular consequence since the contributions of the 24th 
shell are minute. 

Since at low densities the wave- functions are not reliable, 
the energy minimum that appears in the curve of E/N vs. p (Table 3) 
is not to be trusted. One possible source of error could 


be attributed to the HF; we therefore moved a quite arbitrarily 



TABLE. 4 


in. 

\C 


LOCATION OF THE PEAK OF THE WAVE FUNCTIONS (in fermis) 



6. 69 x 

10 12 

10. 

5. 7 

6. 6 

2. 43 x 

10 13 

6. 5 

3. 5 

4.48 

5. 35 x 

10 13 

5. 

1. 85 

3.45 

h - 4 

• 

O 

X 

10 14 

' 4. 

2. 28 

2. 88 

2.84 x 

10 14 

2. 866 

1. 805 

2. 15 

3.96 x 

10 14 

2. 566 

1. 64 

1. 95 


7. 5 

7. 07 

1. 5 

1. 2 

5. 7 

4. 87 

4. 59 

1. 56 

1. 36 

4. 09 

3.75 

3. 53 

1. 65 

1. 35 

3. 3 

3. 12 

2. 82 

1. 68 

1. 32 

2. 76 

2. 23 

2. 02 

1.46 

1. 20 

2. 06 

1. 95 

1. 815 

1.44 

1. 13 

1.85 


5. 1 

4.8 

1. 2 

1. 5 

6. 6 

3. 9 

3. 7 

4. 48 

1. 56 

4.48 

3. 15 

2. 85 

3.45 

1. 5 

3. 45 

2. 52 

2.4 

2. 76 

1.44 

2. 88 

1. 89 

1.80 

1. 97 

1. 29 

2.06 

1. 745 

1.64 

1. 85 

1. 23 

1. 85 
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into two opposite directions: in both cases the results either did 

not change appreciably or became more disastrous. A close 
inspection of Table 4 shows that the misbehavior in the 8- equations 


is much more pronounced than in the system of 3. Among the 
3 3 3 3 

triplets waves, P^, P^» P^ and still behave reasonably 

3 3~ 3 3 ~ 

well whereas S^, S^, and grossly misbehave. The 

reason can be seen to be the uneven presence of the tensor force 
3 3 3 3 

in Eq. (12): P^, P^» P^ and do not contain in the 


same way as the others do; we therefore arbitrarily put = 0 
at low density. The situation improved quite remarkably. The 


h' s were having the maximum at the expected positions. This 

procedure seemed to indicate that the tensor force should contain 

a density dependent factor, say f (p ) , equal to one at p = o 

r nucl. 

which decreases when p decreases. Since only one point on 

the curve is known, p = p . , one cannot hope to determine 

nucl. 


the form of the function. The fact that at low density the values j 

of E/N are misbehaving is not at all surprising if we consider 

i 


that the explicit form of the nucleon- nucleon potential used here 
has been fitted to data at the nuclear density and not any lower. 
The argument is of general nature and one therefore would be 


led to think that not only V T is affected, but all the remaining 
parts of the potential. This is surely true, but a low densities 
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the s-waves are dominant and their misbehavior is more important 

than any other wave; consequently the tensor force that couples 
3 3 

and plays a particularly significant role. 

HCP Structure 


A more interesting case is presented by the HCP structure. 
It is known from solid state physics that in many instances when 
an FCC structure ceases to exist as a result of some of its elastic 


constants going soft, the crystal structure most frequently does 
not melt but, rather, shears into an HCP structure. It was there- 
fore natural to consider the possible alternative. The geometry 
is shown in Fig. 19. The energy results are given in Table 5. 


TABLE 5 


! 


p 

r HCP 

A 

-1 

a 

K. E. 

P. E. 

E/N 

E/N (V T = 

(gr/cc) 

(f) 

(f) 

(£) 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

3.1 X 10 13 

4. 243 

4. 243 

1. 57 

12. 56 

-12.06 

+ 0.4998 

8.48 

13 

5. 35 x 10 

3. 536 

3. 536 

1. 334 

17.43 

-17. 51 

-0. 0811 

10. 60 


1. 04 

X 

10 14 

2. 828 

2. 828 

1.067 

27. 23 

-27. 35 

-0. 1232 

2. 84 

X 

10 14 

2. 026 

2. 026 

0. 765 

53. 04 

-54. 55 

-1. 5 

3. 96 

X 

10 14 

1. 8151 

1. 8151 

0. 672 

68. 64 

- 66. 42 

+ 2. 219 


TABLE 6 


LOCATION OF THE PEAK OF THE WAVE FUNCTIONS (in fermis) 



p 

(gr/cc) 

r HCP 

(f) 

ls o 

lp ! 


A 

(f) 

3 S- 

3 s 

3 -0 

3p i 

*», 


3 5, 

3D 2 

6. 69 y 

10 12 

7. 07 

5. 51 

6.78 

7. 42 

7. 07 

1. 7 

1. 27 

5. 73 

5. 3 

5. 1 

6. 6 

1. 5 

6. 6 

3.1 x 

10 13 

4. 24 

3. 4 

4. 1 

4.4 

4. 24 

1. 7 

1.4 

3. 9 

3. 6 

3.4 

4. 1 

1. 5 

4. 1 

5.35 X 

10 13 

3, 53 

2.7 

3. 5 

3. 8 

3. 53 

1.7 

1. 3 

3. 3 

3. 1 

2.9 

3. 5 

1. 5 

3. 5 

1 . 04 X 

10 14 

_ 14 

' 2. 83 

2. 26 

2. 83 

3. 05 

2. 83 

1. 58 

1. 24 

2.7 

2. 6 

2.4 

2.7 

1. 5 

2. 83 

2. 84 x 

10 

2. 03 

1. 78 

2. 1 

2. 3 

2. 03 

1.46 

1. 22 

2. 03 

1.9 

1. 8 

2.03 

1. 3 

2. 1 

3. 96 x 

10 15 

1. 815 

1.66 

1. 89 

2. 03 

1, 81 

1.45 

1. 16 

1. 89 

1. 74 

1. .6 

1. 815 

1. 3 

1. 89 
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Again the situation becomes unrealistic as soon as the h r s 
start peaking at an unphysical distance. We have also included the case 
= 0. The situation regarding the location is markedly improved 
at very low density, whereas it is probably over-estimated in the 
intermediate region where the join-up should be smoother. In Table 
6 we report the location of the peaks for the various h's. The HCP, 
undoubtedly presents a distinctive characteristic with respect to 
BCC and FCC. For the density region where the wave-functions are 
still acceptable, namely in the lower part of the table, the energy 
E/N vs. p has a minimum at exactly the nuclear density. It is 
very slightly bound -1. 5 MeV. 

IV. ELASTIC CONSTANTS ■ ; 

We have stressed in the Introduction that the system of 
neutrons and protons arranged in a lattice structure must be 
tested for stability against deformations. There have been 
several attempts to derive criteria for melting of a solid, the 
well known amongst which is Lindemann's rule. This rule is 
based on the assumption that a solid melts when the amplitude 
of oscillation of a vibrating particle is a sizable fraction of the 
nearest neighbor distance. Such an empirical criterion can only 
be regarded as a convenient ’’rule of thumb" designed to test 
semi-quantitatively the stability of a crystalline structure and 
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preceding page blank, not filmed 


indeed, a complete theory of melting must necessarily examine 
the detailed stability of a lattice when it is deformed under 
shearing stresses. 


Any lattice has to satisfy the requirement that its energy 
density must have a stationary value at equilibrium. However, 
for the structure to be stable, the energy must have a positive 
definite quadratic form, thus leading to a$ increase in its value 
when undergoing a small strain. For the purpose of investigating 
the mechanical stability of a crystal lattice we consider a cube 
with side 2a which is deformed homogeneously into a parallelepiped 
whose sides are the vectors r , r, , r . Thus the lattice points 

i ^ J 

of the undistorted cube which are described by 


It 

o 


(b a> 


V 


V 


are now given by the vectors 


R 6 " t l r i + l Z X 2 + S r 3 


where the integers assume different sets of values for 


the B. C. C. , F. C. C. lattices. Now for small deformation the 


square of the distance from the origin is given by 
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2 2 2 
R = R + a 
o 


£ 

as 


as 


4 

a p 


n 2 

R + 5 « 
o 


( 16 ) 


where e are the strain components of the theory of elasticity: 

exp 






etc. 


The total potential energy E of a lattice where the interparticle 
forces are only central, may be written as 

E/N = -iE® (Tr^Ts ) 


But since for any function §(r) we have 



and hence the potential energy can be written as 
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(E/N) fi = | * + iD? + 6 2 D 2 i + . . . j 

= (E/N) + f~ y~l e I l Dg + TZ 

o 4 a--/ afi a fi 16 ' 


(17) 


e e £ £ £ £ D $ + ♦ » • 

ag pv a 3 p v 


Here the first term on the right-hand side is the energy of the 
undeformed cube which will be neglected hereafter; the second 
term is zero at equilibrium while the last term may be rewritten 


as 


(E/NV) = ~ c ..(e 2 + e 2 + e 2 \ -f c L 
§ 2 11\ xx yy zz) 12 l 


2 C 44 (' 


2 2 2 

e + e -f e 
yz zx xy 


) 


with 


11 


^ Ed 2 U, 4 -p, 

V 1 


ee + e e + e e 
yy zz zz xx xx yy i 


(18) 


2a ^ 

Y 1 


2 . 2.-2 


12 v 


. 2a Vn 2 , 2. 2 _ 

c 44 T ^ * K \ K Z ' 


y = 4(BCC) , y = 2(FCC) , V = ya , 


and the summation extends over successive shells. 

If we now identify R with A^ » the distance from the origin 
of various shells, we can compute, numerically, the first and 
second derivatives of the potential energy at successive shell 
distances and sum over a desired number of shells (in order to 
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achieve convergence) and evaluate the elastic constants c^ and 

C 12 aS a ^ unc ^ on ma tt e ^ density p. The crystalline 

structure is considered stable if (E/NV) > 0, i. e. , if there is 

6 

a gain in energy while undergoing a deformation. It is straight- 
forward to show that the quadratic form (18) is positive- definite 
provided the following conditions are satisfied: 


C ll + 2C 12 >0 * 


44 


> 0 , 


^11 " C 12 > ° * 


Physically it is clear that the elastic constant c^ represents 

the shear modulus as it multiplies that part of the deformation 

energy which is proportional to the non- diagonal terms, e , 

yz 


zx 


For computing the elastic constants we follow the procedure 

( 8 ) 

outlined by Born 1 ; and the results for nuclear matter arranged 

35 2 

in an FCC lattice are (in units of 10 dynes/cm ): 


11 

-2.83 

C 44= - 71 

r FCC ~ 2 * 57f ' 

p = 3. 96 x 

11 = 

-1. 79 

c 44 = -48 

r FCC = 2f 86f 1 

p = 2. 84 x 


14 

14 


cm 


thus indicating a complete mechanical instability of the system. Here the 


pressure P is very small compared to the C 


o and hence C, „ ss C 

0,0 12 44 . 
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For a BCC structure the corresponding elastic constants 
again come out to be negative and the system is evidently mechani- 
cally unstable. We thus conclude that nuclear matter is unlikely to 
fexist in any of the commonly occurring crystalline structures. . 

V. CONCLUSIONS 

None of the most common solid-state structures seem to be 

an even probable candidate to represent what is known as nuclear 

matter. The conventional BCC and FCC configurations very likely 

do not even give a minimum in the energy curve, or, if they do so, 

it would be in a density region much lower than the usual nuclear 

matter density. However, HCP does have a minimum at 

14 3 

p = P nuc ^ = 2. 842 x 10 gm/cm with a modest binding of energy of 
- 1. 5 MeV, about a factor of ten lower than the experimental value. 

In conclusion, a solid-state structure does not seem to be 
a better alternative to the usually employed gas computation. 

It remains to be seen what an actual liquid type of computa- 
tion would do. Due to the great difficulties already encountered 
in the theory of classical liquids, an extension of that formalism 
to include complicated angular and spin dependent nuclear forces 
does not seem to be just around the corner. 
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FIGURE CAPTIONS 


Fig. 1: 
Fig. 2: 
Fig. 3: 

Fig. 4: 
Fig. 5: 
Fig. 6: 
Fig. 7: 

Fig. 8: 
Fig. 9: 
Fig. 10: 


The B.C.C. arrangement for an equal number of 
neutrons and protons 

The F. C. C. arrangement for an equal number of 
neutrons and protons 

The singlet -6=0 wave -functions h vs. r at a give 
density. 

The peak is very pronounced and it is around the 
first neighbor distance. 

Same as in Fig. 3 for the triplet, -t = 0 case.' 

Same as in Fig, 3 for the triplet, t = 1 case. 

Same as in Fig. 3 for the triplet, t - 2 case. 

The product hV vs. r [ A (*S ) = h etc. 1 for 

.■ f . Y o o 

singlet even. 

The product hV vs. r for the mixed terms. 

The product hV vs. r for the triplet odd case. 

The product hV vs. r for the triplet even case. 


Fig. 11 - 18: The same as for Fig. 3-10 for the 24th shell. 

It can be seen that the wave- functions lag behind, 
i. e. , they do not peak quite at the appropriate A • 
This brings too much attraction into the energy. 
Fortunately the contribution of such a distant shell 
is negligible. 


Fig. 19: 


The H.C.P. configuration for equal number of 
neutrons and protons. 
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